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Abstract 

We present a study of the Dirac theory of constraints emphasizing the du- 
aUty between the Poisson-algebraic and the geometric points of view, related 
respectively to the work of Dirac and of Gotay and Nester. The global notion 
of second class submanifold as being tangent to a second class subbundle is 
studied and equations of motion are written. 

Jerry Marsden passed away on September 21, 2010. He was a friend and a great 
man and mathematician. It was a privilege that he was a cohaborator in this paper. 



1 Introduction 

In this brief introduction we give a sketch of the present work. References, prehmi- 
naries and a description of the main works related to the present paper will be given 
in section 2. 

^Corresponding author. Fax: +54-221-424-5875 
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The Dirac theory of constrained Hamiltonian systems was written by Dirac in 
terms of the canonical Poisson brackets in the space of classical observables (func- 
tions on the phase space) avoiding the notion of constraint submanifolds, which, on 
the other hand, is naturally present in the theory. For instance, instead of intro- 
ducing the final constraint submanifold the notion of weak equality of functions is 
preferred. This indicates a deliberate decision to study the Poisson-algebraic aspect 
of constrained Hamiltonian systems, and for good reasons, in view of the success of 
the theory. 

However, from the point of view of classical mechanics the states are points in 
phase space and the constraint submanifolds also play an important role in under- 
standing the geometry of equations of motion and solution curves. This point of 
view also suggests, for instance, that one can view a constrained Hamiltonian system 
as an IDE (Implicit Differential Equation). Then the Dirac algorithm, as well as 
several questions about dynamics like existence of solutions for the initial condition 
problem, has a meaning also in the context of IDEs, which is an active field of study. 

In many papers after Dirac's work, for instance the work of Gotay and Nester, 
cited in the next section, the geometric side of his theory has been developed and 
proved very useful. The geometric side is, in a sense, dual to the algebraic side 
and this duality is apparent in the commutative diagrams of subsection 3.2 where 
geometric diagrams have a Poisson-algebraic counterpart. This algebro-geometric 
possibility of approaching questions is present throughout the paper. Even though 
there are many beautiful works emphasizing the Poisson-algebraic or the geometric 
aspects our main references will be the works of Dirac and those of Gotay and 
Nester, respectively. 

The notion of Dirac structure (Courant and Weinstein [1988]; Courant [1990]; 
Bursztyn and Crainic [2009]) is originated in part in Dirac's work and gives a new 
possibility to understand and extend the theory. The present paper should be fol- 
lowed soon by a generalization in the realm of Dirac geometry. 

In section 2 we review the Gotay-Nester and the Dirac algorithms and define 
the notion of secondary constraint submanifolds, in particular, the final constraint 
submanifold. Both the primary and the final constraint submanifolds are important 
for writing equations of motion. 

In subsection 3.2 we perform a careful study of the primary and final con- 
straint submanifolds and various quotient manifolds and commutative diagrams, 
which helps to understand important aspects of the dynamics, for instance, the 
question of uniqueness of solutions. The latter is related to the notion of physical 
variables. We also study the dual point of view using the Poisson algebra of first 
class constraints and various quotients and commutative diagrams, which shows the 
duality between the geometric and the Poisson-algebraic points of view. In subsec- 
tion 3.3 we show that the second class constraint submanifolds are submanifolds 
of the phase space tangent to a second class vector subbundle along the final con- 
straint submanifold, which may lead to a classification, at least in some examples, 
of second class constraint submanifolds modulo tangency. The second class vector 
subbundle carries enough information to write the Dirac bracket at points of the 
final constraint submanifold. We write equations of motion in terms of the Dirac 
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bracket in subsection 3.4. 

Another feature of our work is that notions such as second class constraint 
submanifolds and Dirac brackets are defined globally (in the sense explained in 
theorem 3.20), and notions such as the standard Dirac bracket appear as coordinate 
expressions of a global object. 

Along this paper we assume strong regularity conditions that lead for instance, 
to the fact that the Dirac bracket is locally constant, which can be established using 
the Weinstein splitting theorem. 

2 Constraint Algorithms. 

Implicit Differential Equations We now briefly review some basic results con- 
cerning general IDEs and constraint algorithms. Let M be a given differentiable 
manifold. An IDE on M, written as 

¥p(x,x)=0, (2.1) 

of which ODEs x — f{x) = or algebraic equations (t){x) = are considered trivial 
particular cases, appear naturally in science and technology. A solution of (2.1) 
at a point x is a vector {x,x) G T^M satisfying (2.1). A solution curve, say x(t), 
t € {o,,b), must satisfy, by definition, that {x{t),x{t)) is a solution at x{t) for all 
t G {a,b). In the local case, M is an open set of M" and the IDE can be written 
equivalently in the form 

X = u 

= ip{x, u), 

that is, 

A{y)-y = f{y), 

where y = {x,u) and f{y) = {u,f{x,u)). This is by definition a quasilinear IDE. 

Basic questions such as existence, uniqueness or extension of solutions are not 
completely answered yet, although many partial results in this direction have been 
established. One of the common features of those results is that they show, at least 
under certain regularity conditions, how to transform, using a certain constraint 
algorithm, a given IDE into an equivalent parameter dependent ODE on a certain 
final constraint manifold. 

Very briefly, the idea consists in finding a decreasing sequence of constraint 
submanifolds M ^ Mi 5 ... 5 Mc, which appears naturally by imposing the 
condition of existence of a solution (x, x) € TM^ to the given IDE at a each point 
X € M/c_|_i. Under certain additional conditions, like locally constant rank conditions, 
the original IDE is reduced to an equivalent ODE depending on parameters on 
the final constraint submanifold Mc, which, by construction, has the fundamental 
property that it must contain all solutions curves of the given IDE. 

In spite of the simplicity of the general algorithm, there are fundamental exam- 
ples where extra meaningful structures are used to build the submanifolds M^ and 
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to write the equations of motion. Moreover, in the Dirac approach the sequence of 
submanifolds is not emphasized and the methods of Poisson geometry are used, 
for good reasons. 

We can compare Dirac [1950]; Gotay, Nester, and Hinds [1978]; Rabier and 
Rheinboldt [1994]; Pritchard [2003]; Cendra and Etchechoury [2006], to see how the 
idea of the algorithm works in different contexts. In Cendra and Etchechoury [2006], 
one works in the realm of subanalytic sets; in Gotay, Nester, and Hinds [1978] one 
works with presymplectic manifolds; in Pritchard [2003] one works with complex 
algebraic manifolds; Dirac [1950] uses Poisson brackets; in Rabier and Rheinboldt 
[1994] some degree of differentiability of the basic data is assumed, and, besides, 
some constant rank hypothesis is added, essentially to ensure applicability of some 
constant rank theorem. Some relevant references for general IDEs connected to 
physics or control theory, which show a diversity of geometric or analytic methods 
or a combination of both are Carihena and Ranada [1984]; Gracia and Pons [1991, 
1992]; Mendella, Marmo, and Tulczyjew [1995]; de Leon and de Diego [1996]; Ibort, 
de Leon, Marmo, and Martin de Diego [1996]; Marmo, Mendella, and Tulczyjew 
[1997]; Delgado-Tellez and Ibort [2003]. 

In the present paper we will concentrate on the Dirac and the Gotay-Nester 
points of view (which represent the algebraic and the geometric side), see Dirac 
[1950, 1958, 1964] and Gotay, Nester, and Hinds [1978]. One may say that some 
aspects of Dirac's idea have been nicely formalized and generalized in Gotay, Nester, 
and Hinds [1978] in the context of presymplectic geometry on reflexive Banach man- 
ifolds. Both the algebraic and the geometric aspects of Dirac's theory have been 
treated by many people, cited below, with different ideas. The Dirac algorithm is 
not the same as the Gotay-Nester algorithm although the two methods are essen- 
tially equivalent in fundamental examples, like degenerate Lagrangian systems, as 
shown in Gotay, Nester, and Hinds [1978]. The Dirac algorithm provides explicit 
equations of motion written in terms of the canonical bracket of the ambient sym- 
plectic manifold and a total Hamiltonian depending on parameters. Besides, the 
Dirac approach yields the Dirac bracket defined in a neighborhood of the final con- 
straint submanifold. Equations of motion written in terms of the Dirac bracket are 
specially simple and elegant, as we will see. On the other hand, the IDE obtained on 
the final constraints submanifold by the Gotay-Nester algorithm does not depend 
on any parameters or an embedding in a symplectic manifold. 

General assumptions. All manifolds involved will be finite-dimensional smooth 
manifolds and all maps will be smooth, unless otherwise specified. Several argu- 
ments in this paper are of a local character, but they can be regarded as coordinate 
versions of global results. For instance, this is the case for the notion of second class 
constraints, which represent a second class constraint submanifold. 

The Gotay-Nester and the Dirac constraint algorithms considered in this paper 
can be considered as particular ways of writing the general constraint algorithm for 
(quasilinear) IDEs mentioned above, using the special structure available in each 
case (presymplectic structures, symplectic and Poisson structures, respectively). 
Therefore, the sequence of secondary constraints M^, k = 1, . . . is the same for 
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all these algorithms. It is important to have a criterion to ensure that this sequence 
stops. We will assume that each Mk+i is a closed submanifold of Mk defined by equa- 
tions. Also, we assume that for each k and each x G M^+i, if dim^; Mk = diuix M^+i, 
then X G M^+p and dim^; = dima; M^+p for all p € N. This implies that the se- 
quence stabilizes, that is there is some c such that = Mc+p, for all p € N. 
For each one has the corresponding k-th. IDE, which can be written as 

Ak{x) - x = fk{x), 

where Ak{x) = A(x)\TxMk and fk = f\Mk- A point x will be in M^+i, by definition, 
iff this equation has a solution {x,x). We will assume throughout the paper that the 
rank of A^^x) is locally constant on M^+i- This implies that the rank of Ac{x) is 
locally constant on Mc, and on each point of Mc there is at least one solution that is 
tangent to Mc- The main property of Mc from the dynamical point of view is that 
every solution curve to the original system must lie on Mc- Since we are assuming 
the locally constant rank condition, the final system on Mc can be converted, at 
least locally, into a parameter-dependent family of ODEs. 

2.1 A brief review of Dirac's theory. 

Dirac's theory of constraints has been extensively studied from many different 
points of view and extended in many directions. An important part of those 
developments is contained in Sudarshan and Mukunda [1974]; Gotay and Nester 
[1979]; Mukunda [1980]; Gotay and Nester [1980]; Skinner [1983]; Skinner and Rusk 
[1983a, b]; Carihena and Rahada [1984]; Carihena, Gomis, Ibort, and Roman [1985]; 
Cantrijn, Carihena, Crampin, and Ibort [1986]; van der Schaft [1987]; Mukunda 
[1987]; Carihena, Lopez, and Roman-Roy [1988]; Carihena [1990]; Henneaux and 
Teitelboim [1992]; Krupkova [1994]; de Leon, de Diego, and Pitanga [1995]; Krup- 
kova [1996]; Marmo, Mendella, and Tulczyjew [1997]; Ibort, de Leon, Marrero, and 
Martin de Diego [1999]. 

There is a certain duality between the Dirac approach, in which the role of con- 
straints as being functions on the phase space and the canonical bracket is essential, 
and the approach of many other authors, where, in addition, the geometry behind 
the canonical Poisson algebra on phase space is emphasized. This duality is present 
along this paper, and in this sense, our main references will be Dirac [1964] and 
Gotay, Nester, and Hinds [1978]. 

We will recall some essential aspects of the Dirac theory of constraints, follow- 
ing Dirac [1964], but using a more modern language, adapted to our purposes, and 
assuming explicitly certain regularity conditions. 

Dirac's theory starts with a given singular Lagrangian system L : TQ — )• M, since 
in the case of a regular Lagrangian the theory becomes trivial. Then, in order to 
construct a Hamiltonian theory one must consider from the beginning the image of 
the Legendre transformation, which may be a very complicated subset Mq C T*Q. 
By definition, solution curves to the Hamiltonian system must be exactly the image 
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of solutions to Euler-Lagrange equations under the Legendre transformation. Nat- 
ural questions such as existence of solution curves to the Euler-Lagrange equations 
for a given initial condition are not completely solved in general, to the best of our 
knowledge. In order to obtain this kind of results one would need to choose mathe- 
matically precise hypotheses, a topic not considered in Dirac's work. 

In this paper we will assume a general hypothesis about regularity, under 
which this kind of problem is easier. Regularity means, among other hypotheses to 
be established along the paper as they are needed, that certain sets Mq I) Mi . . . 
are submanifolds of T*Q defined regularly by equations (pf^^ = 0, i = l,...,afc, 

(k) 

k = 0, 1, ... . The (/)• are functions defined on T*Q constructed by the Dirac 
constraint algorithm and called constraints. The submanifolds are called 
the constraint submanifolds. 

As usual, one assumes that the algorithm stops for k = c. One also assumes that 
the ranks of the matrices {{(pf^^^ (l>^^^}{x)) and {{4'i^\4'f''^}{^)) locally constant 
on the final constraint submanifold Mc. 

Each (j)f''^ is called a primary constraint and each (p^'^^ is called a final con- 

(k) 

straint. Generically, (j)- , k = l,...,c, are called secondary constraints. The 
main property of the final constraint submanifold Mc is that any motion of the 
classical particle, that is, any solution {q{t),p{t)), must remain in Mc, and Dirac 
shows how to write Poisson equations of motion in terms of position and momentum 
using the canonical Poisson bracket on T*Q and the total Hamiltonian Ht- For 
given initial conditions belonging to Mc, solutions are not necessarily unique and 
Dirac interprets this fact as being due to the nonphysical character of some of 
the variables. Of fundamental importance for Dirac's theory, especially for quan- 
tization, are the classification of constraints into first class and second class in 
terms of certain commutation relations, and the construction of a Poisson bracket 
{ , }* called Dirac's bracket. An important result is that with respect to the Dirac 

(c) 

bracket all final constraints appear to be first class constraints, in other words, 

{cl)[^\ (p^^^ }* {x) = 0, for all x € Mc. Dirac's procedure also shows how to deal with 
the nonphysical variables and find the correct notion of state of the system. One 
shows that there are physically meaningful variables in terms of which the evo- 
lution for a given initial state is determined. This is important from the classical 
and also from the quantum mechanics point of view. 

Now, we will be more precise. The image of the Legendre transformation FL : 
TQ — > T*Q, that is, Mq = FL(TQ), contained in the canonical symplectic manifold 
T*Q, is assumed to be defined by equations (/)-°^ = 0, i = 1,... ,ao, where each 
^(0) : _). ]^ is a primary constraint, by definition. 

In the case in which the Legendre transformation is degenerate the Hamiltonian 
H : T*Q — )• M is not uniquely defined from the formulas pv — L{q,v) = H{q,p), 
p = dL(q,v)/dv, but in Dirac's theory one assumes that such a function can be 
conveniently defined on Mq (which can be done in examples using ideas akin to 
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the Pontryagin maximum principle, like the fact that for each {q,p) G Mq the 
derivative of pv — L{q, v) with respect to v is 0) and then extended, more or less 
arbitrarily, to T*Q. Then one defines, following Dirac, the total Hamiltonian Ht = 
H + X'f^Q-^(f)f\ with arbitrary parameters A^q^ to be determined. The Dirac constraint 
algorithm goes as follows. The preservation of the primary constraints is written 
{cp,f\HT}{x) = 0, i = 1, . . . , ao, X G Mo, or 

{(pf\H}{x) + \\^^{(pf\(t>f}{x) = 0, i,j = 1, . . . ,ao, X G Mo. 

Then Mi is defined by the condition that x G Mi if and only if there exist A(o) = 

(A^Q^, . . . , A"q^) such that the system of equations 4)f'\x) = 0, {(l>f'\ Ht}{x) = 0, i = 
1, . . . , ao) is satisfied. Clearly, Mo 5 Mi, and one assumes that Mi is a submanifold 
regularly defined by equations, say, (j)^^'^ = 0, i = l,...,oi, where each c/)^^-' is a 
secondary constraint, by definition. By proceeding iteratively one obtains a sequence 
Mq 2 Ml 2 • • • , and we will assume that this sequence stops. Then there are final 
constraints, say i = 1, . . . ,ac, defining regularly a (nonempty by assumption) 

(c) 

submanifold Mc by equations 4>l = 0, i = l,...,ac, called the final constraint 
submanifold, and the following condition is satisfied. For each x G Mc there exists 
(A^Qj , ■ ■ ■ , ^lo)) such that 

{<P^\h}{x) + Ajo){4'\</.f }(x) = 0,i = l,...,a„j = l,...,ao. (2.2) 

For each x G Mc the space of solutions of the linear system of equations (2.2) 
in the unknowns A|q^ is an affine subspace of M'^o, called si^^ whose dimension is 

a locally constant function d^'^\x) = oq — Taiak{{(j)\^\(p^^^}{x)). One can locally 
choose d^^\x) unknowns as being free parameters and the rest will depend affinely 
on them. Then the solutions of (2.2) form an affine bundle S^'^^ over Mc. After 
replacing A(o) G S^^^ in the expression of the total Hamiltonian, the corresponding 
Hamiltonian vector field, 

Xht{x) = Xh{x) + A^g)X_^(o)(a;), 

X G Mc, which will depend on the free unknowns, will be tangent to Mc. Its integral 
curves, for an arbitrary choice of a time dependence of the free unknowns, will 
be solutions of the equations of motion, which is the main property of the final 
constraint submanifold Mc from the point of view of classical mechanics. The lack 
of uniqueness of solution for a given initial condition in Mc, given by the presence 
of free parameters, indicates, according to Dirac, the nonphysical character of some 
of the variables. In our context the physical variables can be defined on a quotient 
manifold. 

Remark Dirac introduces the notion of weak equality for functions on r*Q. Two 
such functions are weakly equal, denoted f ^ g, if /|Mc = g\Mc. Then, for instance 

(k) ■ (c) 

RiO. If / then / = i^'cj)] ', for some functions z^* on T*Q and conversely. 
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Since we have introduced the notion of a constraint submanifold, in particular the 
final constraint submanifold, we prefer not to use the notation ~. 

Now let us make some comments on the notions of first class and second class 
constraints. The rank of the matrix {{cl)f^\(j)^^^^}{x)),i,j = 1, . . . ,ac, is necessarily 
even, say, 2s, and it is assumed to be constant. Then, using elementary properties of 
determinants (like adding to a row or column a linear combination of the other rows 
or columns) one can find, at least locally in a neighborhood of each point x € Mc, 
functions ipi, i = 1, . . . , — 2s, with, Xji j = li • • • j 2s, such that the equations 
ipi = 0, Xj = 0) define Mc regularly and, besides, {ipi,ipi'}{x) = 0, {ipi,Xj}{x) = 0, 
det({xj, Xj'}(3^)) 7^ 0) for hi' = ^, ■ ■ ■ ,cic — 2s, = 1, ... ,2s and x G M^.. In 

(c) 

fact, we will assume that this is can be done globally, for simplicity. The <j))j are 
linear combinations with smooth coefficients of the Xj s-^id ipi, and conversely. The 
functions J = !> • • • ) 2s, are called second class constraints and the functions 
■^j, i = — 2s, are called first class constraints. 

More generally, any function p on T*Q satisfying p\Mc = 0, {p,ipi}\Mc = 0, 
{p,Xj}\Mc = 0, is a first class constraint with respect to the submanifold Mc, by 
definition. Any function g on T*Q satisfying {g,tl'i}\Mc = 0, {g,Xj}\Mc = 0, is a 
first class function, by definition. For instance, the total Hamiltonian Ht is a 
first class function. 

Now define the Hamiltonian he in terms of ipi, Xji ^ = 1, . . . ,k, j = 1, ... ,2s, as 

hc = H + X'iPi + p^Xj- 

The preservation of the constraints for the evolution generated by he can be rewritten 
as {ilji,hc}{x) = 0, which is equivalent to {ipi, H}{x) = for all x G Mc, and 
{Xj, hc}{x) = 0, for all x G Mc- The latter is equivalent to 

{Xi,H}{x) + p^Xi,Xj}{x) = 0,i,j = 1, ... ,2s, 

for all X € Mc, which determines the p^ as well-defined functions on Mc. Then the 
solutions {p{x). A) form an affine bundle with base Mc and whose fiber, parametrized 
by the free parameters A, has dimension Uc — 2s. 

Any section {p{x),X{x)) of this bundle determines he as a first class function. 
This means that Xh^{x) G T^Mc, for each x G Mc, and therefore a solution curve of 
Xh^ is contained in Mc provided that the initial condition belongs to Mc- We will 
show in this section that there is a symplectic manifold Mc such that one can pass 
to the quotient Mc Mc and also he passes to the quotient he in such a way that 
solution curves of he become solutions curves of the Hamiltonian he. Moreover, we 
will show that there is a manifold Mc and natural maps Me Me Mc such that 
the Hamiltonian he passes to the quotient to a function he on Mc- The extended 
Hamiltonian defined by Dirac is related to he. One can show that he passes to the 
quotient via Me Me to the function he defined above. 

Dirac defines an interesting bracket, now called the Dirac bracket, 
{F,Gr = {F,G} - {F,x^V'{XJ,G}, 
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which is defined on an open set in T*Q containing Mc, where c*-', which by definition 
is the inverse matrix of {Xi^Xj}i is defined. The Dirac bracket is a Poisson bracket 
and has the important property that, for any function F on T*Q, the condition 
{F, Xj}* = 0, j = 1, . . . , 2s, is satisfied on a neighborhood of Mc, which imphes that 
F = {F, he} = {F, /ic}*, for any function F. Besides, {^j, ^j}* = 0, i, j = 1, . . . , /c, 
on Mc- Because of this, one may say that, with respect to the Dirac bracket, all the 
constraints Xj, J = li • • • j 2s and ipi, i = 1, . . . ,k, are first class with respect to Mc- 
This is important for purposes of quantization. 

2.2 A brief Review of the Gotay-Nester Theory 

In this section we recall some aspects of the Gotay-Nester theory Gotay, Nester, 
and Hinds [1978] which we need. This theory was developed in Gotay, Nester, and 
Hinds [1978] to deal geometrically with the Dirac-Bergman theory of constraints. 
The main equation studied is an IDE of the type 

ii = a(x), (2.3) 

where w is a closed 2-form on a manifold M and a G Q^{M) is a closed 1-form on 
M. As we have indicated before this kind of equation appears naturally in classical 
Lagrangian mechanics, in fact, we will show later that the Euler-Lagrange equations 
can be rewritten equivalently in the form 

ii a;(x) = d£{x), (2.4) 

which is clearly of the type (2.3). 

Description of the Gotay-Nester Algorithm. As we mentioned in the Intro- 
duction, in order to deal with IDEs one can apply a basic idea which consists in 
building a sequence of constraint submanifolds. 

Let us first describe that basic approach for a system like (2.3) without us- 
ing explicitly the presymplectic form, and later on we will briefly explain how the 
presymplectic form can be used to write equations for the constraints submanifolds 
explicitly. The latter is an important contribution of the Gotay-Nester algorithm. 
We will assume for simplicity that M is a vector space. 

We want to find solution curves to (2.3). Let x{t) be such a solution curve; then 
for each t the linear algebraic system 

i„(j) a; (x(t)) = a{x{t)), 

has at least one solution, namely, v{t) = x{t). This implies that, for each t, x{t) 
must belong to the subset 

Ml = {x € M I it, u}{x) = a{x) has at least one solution v G T^M}. 

Assume, as in Gotay, Nester, and Hinds [1978], that Mi is a submanifold of M. 
Since x{t) € Mi for all t we must have that x{t) E T^^^^Mi for all t. This implies 
that, for each t, x{t) must belong to the subset 

M2 = {x £ Ml I iyU}{x) = a{x) has at least one solution v € T^Mi}. 
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We can continue in a similar way and define Mfc_|_i recursively as 

Mfc+i = {x € Mk I \vIjj{x) = a{x) has at least one solution v € TjM^}. 

This sequence stabilizes, under the General Assumptions described in section 1. 
Under the assumption that the map 

J : T^Mc T*M\Mc 

has locally constant rank on the final constraint manifold Mc, existence of local 
solution curves to (2.3) for each initial condition in is guaranteed. For given 
local coordinates (xi, . . . ,Xm) on Mc and for a given initial condition xq € Mc 
one can fix some appropriate coordinates as functions of t, say xi{t), . . . ,Xr{t), 
where r = dimkerw, in a neighborhood of xq and then solve (2.3) uniquely for 
Xr+i{t), . . . ,Xm{t). More precisely, in local coordinates our equation becomes 

ujij{x) ■ x^ = ai{x), (2.5) 

i = 1, . . . ,n, j = 1, . . . ,m. Since r = dimker cj, in a neighborhood of a given point xq 
one can solve, after relabeling the coordinates if necessary, x^+k = Xr+k{xi, ■ ■ ■ , Xr), 
k = 1, . . . ,m — r. After choosing arbitrarily the curves xi{t), . . . , Xr (t) and replacing 
in (2.5) one obtains a time dependent ODE in the remaining variables x^-^i, . . . , x^ji. 
We can also interpret the previous arguments by saying that the implicit differential 
equation 

uj{x){x, ) = a(x), X G Mc 
is an ODE on Mc, depending on m — r parameters. 

The following geometric description will be useful later on. Since by assumption 
r does not depend on x G Mc, at least locally, then the equation on Mc 

u;ix){X,) = a{x), (2.6) 

where X S TMc, defines an affine distribution on Mc of locally constant rank. More 
precisely, one has an affine bundle S^'^'^ with base Mc whose fiber si^^ at a given 
point X G Mc is, by definition, 

5^ = {X £ T^Mc I (2.6) is satisfied}. (2.7) 

Remark, (a) If rank uj{x) is not locally constant we still have a distribution S^^^ 
on Mc, but it may be singular. The analysis of existence of solution curves in 
this case may be problematic, see Cendra and Etchechoury [2006], Pritchard [2003] 
and references therein. The algorithm developed in Cendra and Etchechoury [2006] 
for a general system of the type a{x) ■ x = f{x), with analytic data, represents 
an improvement of the previous basic ideas also in the sense that the final system 
obtained after applying the algorithm has always locally constant rank, and that 
singular cases are also studied using desingularization methods. 

(b) In Gotay, Nester, and Hinds [1978] it is explained how solutions can be 
expressed using brackets, as in Dirac's work. 
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Example. Let L : TQ ^ M be a Lagrangian, degenerate or not. Since the problem 
is of a local nature we can use local coordinates. Let £{q,v,p) = pv — L{q,v) and 
let uj € ^}'^(TQ © T*Q) be the presymplectic form u = dq^ A dpi on the Pontryagin 
bundle M = TQ @T*Q. Then Euler-Lagrange equations are written equivalently 
in the form of equation (2.4) with x = {q,v,p). In fact, we have 

Hq,v,p)dq' A dpi = q'dpi - pidq' (2.8) 

d£ = ^^dq^ + ^dp, + (2.9) 
oq^ opi ov^ 

= - — {q,v)dq^ + v^dp, + i^p^ - —^j dv' (2.10) 
Using equations (2.8)-(2.10) we can easily see that (2.4) is equivalent to 



which is clearly equivalent to the Euler-Lagrange equations. The idea of using the 
Pontryagin bundle to write important equations of physics like Euler-Lagrange or 
Hamilton's equations appears in Livens [1919]; Skinner [1983]; Skinner and Rusk 
[1983a, b]; Cendra, Marsden, Pekarsky, and Ratiu [2003]; Yoshimura and Marsden 
[2006a,b]. 



Describing the Secondary Constraints Using uj Now we shall explain how 
the constraint manifolds Mk defined by the algorithm can be described using the 
presymplectic form u. This is a simple but important idea because, in particular, 
we show how each Mk is defined by equations written in terms of co. Depending on 
the nature of w one may obtain analytic, smooth, linear, etc., equations, which may 
simplify matters in given examples. This idea is also important in the context of 
reflexive Banach manifolds, as remarked in Gotay, Nester, and Hinds [1978]. Besides, 
those equations will be obviously invariant under changes of coordinates preserving 

OJ. 

The idea is the following. The condition defining the subsets M^+i, A; = 0, 1, ... , 
calling Mq = M to uniformize the notation, namely, 

i„w(x) = a{x) has at least one solution v € T^M^, 

can be rewritten equivalently as a{x) € (TxMk)^ . Taking into account that (T^Mk)^ = 
{{TxMk)^)° we can conclude that Mfc_|_i is defined by the equation 



Mfc+i = {xeMk\ {a{x), {T^MkY) = {0}}. 
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3 Main Results 
3.1 Preliminaries 

We will need the following results about linear symplectic geometry which are an 
essential part of many of the arguments in our treatment of Dirac and Gotay-Nester 
theories. This is because under our strong regularity assumptions those theories are, 
to a certain extent, linear. 

Lemma 3.1. Let {E,Q) be a symplectic vector space of dimension 2n, V (1 E a 
given subspace. For a given basis ai, i = 1, . . . , r of V° , let Xi = a\, i = 1, . . . r. 
Then the rank of the matrix [ai{Xj)] is even, say 2s, and Xi, i = 1, . . . ,r form a 
basis of . Moreover, the basis ai, i = 1, . . . ,r can be chosen such that for all 

ai{Xj) = 6ij-s, 1 < « < s 
ai{Xj) = -6i-sj, s + 1 <i <2s 
ai{Xj) = 0, 2s <i<r. 

Proof. Consider the subspace = {V°)'^. By a well-known result there is a basis 
Xi, i = 1, . . . , r of such that for all j = 1, . . . , r 

n{Xi, Xj) = 5ij-s, 1 < i < s 
^l{Xi, Xj) = -6i-s,j, s + l <i <2s 
n{Xi,Xj) = 0,2s<i<r 

then take Oi = X'^. The first part of the lemma is easy to prove using this. ■ 

Lemma 3.2. Let ai, i = 1, . . . ,r be a basis of V° having the properties stated in 
Lemma 3.1. Then Xi, i = 2s + 1, . . . ,r form a basis ofVCi V^. 

Proof. Let X = X'Xi be an arbitrary vector in V^. Now X^Xi G V D iS 
aj{X) = X^{aj{Xi)) = 0, j = l,...,r. Since the first 2s columns of the matrix 
[aj(Xj)] are linearly independent and the rest are zero, we must have A* = 0, for 
1 < i < 2s, and A*, i = 2s + 1, . . . ,r are arbitrary. This means that V H is 
generated by Xi, i = 2s + 1, . . . ,r. ■ 

Corollary 3.3. diuiV D = r - 2s. 

Proof. Immediate from lemma 3.2. ■ 

Let CO be the pullback of to y via the inclusion, lo = j*^. Then (V,u!) is a 
presymplectic space. In what follows, the and " operators are taken with respect 
to J7 unless specified otherwise. 

Lemma 3.4. = VnV^. 



Proof. X G iff uj{X,Y) = 0, G F iff n{X,Y) = 0, VY G This is 
equivalent to X eV CiV^. ■ 
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Lemma 3.5. Let ji, i = 1, . . . r be a given basis of V° and let Yi = jf, i = 1, . . . r . 
Let f3 ^ E* be given. Then the following conditions are equivalent. 

(i) /3(y^) = 0. 

(a) The linear system 

^{Y) + X^-fj{Yi) = (3.1) 
has solution X = (A^, . . . , A''). 
Proof. Let us show that (3.1) has solution (A"^, . . . , A'") iff the system 

P{Xk)+fi'ai{Xk) = (3.2) 

has solution (/i^, . . . , fi^), where k,l = 1, . . . r and a/ is a basis satisfying the condi- 
tions of lemma 3.1. Since Yi, i = 1, . . . r and X^, k = 1, . . . r are both bases of 
there is an invertible matrix [^4^ such that Xk = A^Yi. Let [BI] be the inverse of 
[Al.], so Yi = B\Xi. Assume that (3.1) has solution A-', j = 1, . . .r. We can write 
(3.1) as 

/3(yi) + A^o(yj,yi) = o, i = i,...,r. 

Using this we have that for k = 1, . . . , r 

= liiA'kYi) + XmiYj.A'^Yi) = p{Xu) + X^VL{Yj,Xk) 
= p{Xk) + Xin{B\Xu Xk) = /3(Xfc) + ^l'^{XuXk) 

where = B^-. This means that the system (3.2) has solution. The converse is 
analogous. Using this, lemmas 3.2 and 3.4, and the form of the coefficient matrix 
[Q;;(Xfc)] in lemma 3.1, the proof that (3.2) has solution = {p} , . . . , /i^) iff fi{V^) = 
is easy and is left to the reader. ■ 

Lemma 3.6. Consider the hypotheses in lemma 3.5. Then the solutions to 

ixuj = m (3.3) 

(if any) are precisely X = + X^Yj, where {X^ , . . . , A*") is a solution to (3.1). A 
solution to (3.3) exists if and only if /3{V'^) = 0. If uj is symplectic then (3.1) and 
(3.3) have a unique solution and if, in addition, (3^ € V , then A"^ = 0, . . . , A*^' = 
and (5^ coincides with X = {j3\V)^'-' defined by (3.3). 

Proof. Since Yj, j = l,...r form a basis of we have that (A"*^, . . . , A*^') is a 
solution to (3.1) iff (^ + X^-ij){V^) = iff ^ + A^^j G iff ^« + X% G V. Now, let 
X = + A^Yj, where (A\ . . . , A*^) satisfies (3.1). Then we have X & V as we have 
just seen and we also have 

ixu: = {ix^V = X^\V = {P + X^^j)\V = /3|y, 

since 7j, j = l,...r generate V° . We have proven that X is a solution to (3.3). 
To prove that every solution X to (3.3) can be written as before, we can reverse 
the previous argument. Using this, it is clear that if u is symplectic then (3.1) has 
unique solution, in particular, we have that det {'~lj{Yi)) 7^ 0. If, in addition, 
then X^Yj = X — /S" G y. Since 1^-, j = 1, . . . , r is a basis of , using lemma 3.4 
and the fact that = {0} we get that A^ = for j = 1, . . . , r. ■ 
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Corollary 3.7. Let A = {A | A satisfies (3.1)}. Then dimA = r — 2s = dimkerw. 

Proof, kerw = V'^ , which has dhnension r — 2s from corollary 3.3 and lemma 3.4. 
On the other hand the dimension of the subspace of A satisfying (3.1) is clearly also 
r — 2s, since the coefficient matrix has rank 2s. ■ 

3.2 A Poisson- Algebraic and Geometric Study of the Primary and 
Final Constraint Submanifolds. 

The Dirac algorithm, briefly explained in the previous subsection, can be applied 
to any given constrained Hamiltonian system {P, $7, H, M) where (P, Q) is a 
symplectic manifold, the primary constraint submanifold M is a given submanifold 
of P deflned regularly by an equation i;^ = and H is a, Hamiltonian defined on P. 
This is because the particular cotangent bundle structure of the symplectic manifold 
T*Q is not essentially used in the Dirac algorithm. 

For instance, an interesting variant of the Dirac algorithm for a degenerate La- 
grangian system is the following. 

Consider the canonical symplectic manifold N = T*TQ with the canonical 
symplectic form Q, and let the primary constraint be M = TQ ® T*Q, canoni- 
cally embedded in via the map given in local coordinates {q,v,p,v) of N by 
{q,v,p) I— )• {q,v,p,0), in particular, M is defined regularly by the equation u = 0. 

Remark This embedding of M in N is globally defined, and it is a coisotropic 
embedding, see Appendix A. For a given presymplectic manifold (M, a;) one can 
always find a coisotropic embedding in a symplectic manifold [P, Q), see Gotay and 
Sniatycki [1981/82]. The number pv is a well-defined function on M and can be 
naturally extended to a function on a chart with coordinates {q,v,p,v), but this 
does not defines a function on N consistently. In any case, it can be extended to 
a smooth function on N and any such extension will give the same equations of 
motion. The Dirac theory of constraints is essentially a local theory. However, we 
will see a global version of the notion of Dirac bracket, in a sense, as well as its local 
descriptions. 

Consider the function £ : M — )• R given hj £ = pv — L(q,v). Using the fact 
that £ can be extended naturally on a chart with coordinates {q, v,p, u) and taking 
an appropriate partition of unity we can choose once and for all an extension to a 
smooth function £ on N called the Energy. Then we can apply the Dirac algorithm 
to the constrained Hamiltonian system {N,n,£,M). 

There are some interesting features in this approach, as compared with the one 
described in subsection 2.1, where the symplectic manifold is T*Q and the primary 
constraint is the image of the Legendre transformation. For instance, the primary 
constraint M = TQ (B T*Q is well defined in a natural way as a closed submani- 
fold of the symplectic manifold N. Besides, the comparison with the Gotay-Nester 
approach becomes clear from the beginning and the Euler-Lagrange equations are 
derived quickly as a differential-algebraic equation (DAE). On the other hand, this 
approach may have the disadvantage of introducing the extra variable v, which may 
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lead to longer calculations in some examples. 

We shall start with the constrained Hamiltonian system (A^, 0, f , M), and we will 
work locally, for simplicity. We will call (p^^^y i = 1, . . . ,ro, the primary constraints 
z^* defining Mq = M regularly by equations. 

We will emphasize the Gotay-Nester point of view and we will see how it com- 
bines with the Dirac procedure. 

Accordingly, we shall study the Dirac dynamical system on the manifold M = 
TQ (S)T*Q, already considered in the example at the end of section 2.2 where the 
Dirac structure is associated to a presymplectic form w which is the pullback of 
the symplectic form on iV = T*TQ via the inclusion M Q N. Then M is the 
primary constraint. Then the equation to be solved, according to the Gotay-Nester 
algorithm, is the equation 

ioix){X, ) = dS{x)\T,M, (3.4) 

where X € T^Mc and x € Mc, being the final constraint. Let Uc be the pullback 
of via the inclusion of in N . Since ujc is presymplectic, kerwc is an involutive 
distribution. 

From now on we will assume the following. 

Assumption Ki. The distribution kert^c has constant rank and defines a reg- 
ular foliation Kc, that is, the natural map '■ -^c where Mc = Mc/Kc is a 
submersion. 

Lemma 3.8. The following assertions hold: 

(a) There is a uniquely defined symplectic form uJc on Ale such that p*^jjOc = ojc. 

(b) Let X he a given vector field on M^. Then there is a vector field X on 
that is pKc-'^e.lated to X . 

(c) Let f G J-{Mc). Then there exists a vector field X on Mc such that X is pKc' 
related to Xj, and for any such vector field X the equality uJc{x){X, ) = d[p*j^^f){x) 
holds for all x € Mc. 

(d) Let XxQ € TxqMc- Then one can choose the function f € J-{Mc) and the 
vector field X in (c) in such a way that X{xq) = Xxq. 

Proof, (a) By definition, the leaves of the foliation Kc are connected submanifolds 
of Mc, that is, each p^^(z), z G Mc, is connected. For z € Mc, let x € Mc such that 
PkS^) = ^- For A, B € TzMc, as pK^ is a submersion, there are A, B G Mc such 
that TxPKc^ = ^5 TxPKcB = B. We define u)c{z){A, B) = ujc{x){A, B). To prove 
that this is a good definition observe first that it is a consistent definition for fixed 
X, which is easy to prove, using the fact that kera;c(x) = ]ieTTxPKc- Now choose a 
Darboux chart centered at x, say U x V, such that, in this chart, px^ ■ U xV^U 
and C0c{x^,x'^) = Wc(x^), where lOc{x^,x'^) and Cdc{x^) are independent of {xi,X2). 
This shows that Wc is well defined on the chart. Using this and the fact that one 
can cover the connected submanifold pj}{z) with charts as explained above, one can 
deduce by a simple argument that u)c{z) is well defined. 
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(b) Let 5 be a Riemannian metric on Mc- Then for each x E Mc there is a 
uniquely determined X{x) E T^Mc such that X{x) is orthogonal to keiTxPKc s-nd 
TxPKc^{x) = X{x), for all x E Mc- This defines a vector field X on Mc which is 
PK^-^elated to X. 

(c) Given / and using the result of (b) we see that there is a vector field X on 
Mc that is pi^^-related to Xj. Then, for every x E Mc and every E T^Mc, 

uJc{x){X{x),Yx) = u)c{pkXx)) {Xj{pK,{x)),TxPK,Yx) = df{pKAx)){TxPK,Yx) 
= d{p*Kj){x){Yx). 

(d) One can proceed as in (b) and (c) and choosing / such that [df {pKc{xq)))^ = 
TxqPKc^xo stnd, besides, the metric g such that Xx^ is perpendicular to keiTx^^PKc- 

m 

Definition 3.9. (a) For any subspace A C J-{N) define the distribution Aa_ ^ TN 
by AA{x) = {Xf{x)\feA}. 

(b) The space of first class functions is defined as 

i?W = {/ E T{N) j Xfix) E TxMc, for all x E MJ. 

In other words, R^^^ is the largest subset of F{N) satisfying 

^R(c){x) C TxMc, 

X E Mc. 

Remark, (a) From the point of view of classical mechanics, the constraint sub- 
manifolds M and Mc seem to be at least as important as the functions 4>f'^ and 

(f)f^ (called constraints by Dirac) defining them by equations = and (f)f^ = 0, 
respectively. Dirac was interested in classical mechanics, where states are points in 
phase space, as well as in quantum mechanics where functions are observables and 
states are not points in phase space. In the present paper we focus mainly in classi- 
cal mechanics, and therefore we need to concentrate on the constraint submanifolds. 
In particular, M and Mc are the only ones that play an important role. The other 
secondary constraints submanifolds seem to be less important. 

(b) The total Hamiltonian Ht is a first class function, by construction. 

Lemma 3.10. (a) R^^^ is a Poisson subalgebra of {J-{N),{ , }). 

(b) Mc is an integral suhmanifold of Aj^(c). Moreover, for any vector field X on 
Mc that is pK^-related to a vector field Xj on Mc there exists a function f E i?*-'^-' 
such that f\Mc = Px^f '^^d X = Xj. In particular, any vector field X on Mc 
satisfying X{x) E kerujc{x) for all x E Mc is px^-felated to the vector field on 
the symplectic manifold Mc, which is associated to the function / = 0, therefore 
there exists a function f E R^'^\ which satisfies f\Mc = 0, such that X{x) = Xf{x), 
X E Mc. 
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Proof, (a) Let f,g e R^^l Then Xf{x) and Xg{x) are both tangent to Mc at 
points X of Mc which imphes that —X^f gy{x) = [Xf,Xg]{x) is also tangent to Mc 
at points of x of Mc- This shows that {/, g} G R^'^\ It is easy to see that any hnear 
combination of f,g and also fg belong to R^^\ 

(b) By definition Aj^(c) C TMc- We need to show the converse inclusion. Let 
^xo ^ T^qMc, we need to find / € R^^^ such that Xj(xo) = X^^. Choose the function 
/ and the vector field X on Mc as in lemma 3.8, (d). Choose any extension of p*x^f 
to a function g on N . For each x € Mc, we can apply lemmas 3.5 and 3.6 with 
E := T^N, V := T^Mc, /? := dg{x), ji{x) := i = l,...,rc. We obtain 

that in a neighborhood [/ C of each point xq of Mc we can choose C°° functions 
'^(c)(^); ^ = ■ ■ ■ ) '"c; such that 

X(x) = (d<7(3:) + Aj,)(x)d</.f (x))", 

for all X € Mc fl U. Let fu{x) = g{x) + A*^^^(x)0''^^(x), for all x £ U. Then we have 
that X/y(x) = for all x G McCiU. 

Now, consider a partition of unity pi, i € I, on N, where each pi is defined on 
an open set Ui, i E I. Let J C / be defined by the condition i G J if and only if 
Ui n Mc 7^ 0. Using standard techniques of partitions of unity and the above result 
one can assume without loss of generality that for each i € J there is a function fjj. 
defined on Ui such that Xj^ (x) = X{x), for all x G Mc n Ui. Let / = PifUi, 
which can be naturally extended by on A^. Then it is easy to see, using the fact 
that fui{x) = P*K^f{^) = 9(^)5 for each x e Mc, that X{x) = Xf{x), for each 
X € Mc, and in particular Xq = X{xq) = Xf{xQ). ■ 

Lemma 3.11. (a) Each function f E R^^^ is locally constant on the leaves of Kc 
therefore, since they are connected, for each f G R^'^^ there is a uniquely determined 
f S T{Mc), called {pKc)*fj such that f\Mc = Px^f- Moreover, the vector fields 
Xf{x), X G Mc, on Mc and Xj on Mc are pKc-i^slated. 

(b) For each f € J-{Mc) there exists f € R^'^^ such that f\Mc = P*k^I ^.''^d the 
vector fields Xf{x), x E Mc, on Mc and Xj on Mc are px^-related. 

Proof, (a) Let / G R^^\ we only need to show that / is constant on the leaves of 
Kc^ which is equivalent to showing that df[x) \ 'kei ujc{x) = 0, for all x G Mc- For a 
given X G Mc, let G kera;c(x); then using lemma 3.10 (b) one sees that there is 
a function g G i?^"^^ such that Vx = Xg{x). Then we have 

= ujc{x) {Xf{x),Xg{x)) = n{x) {Xf{x),Xg{x)) = df{x)Xg{x). 

Now we shall prove that Xf and Xj are pi^^-related. For each x G Mc and each 
Yx G TxMc, we have 

Uc{Xf{x),Yx) = n{Xf{x),Yx) = df{x){Yx) = d{p*Kj){x){Yx). 

Using this we obtain 

UJc{x){Xf{x),Yx) = Ojc{x){pK,{x)) {TxPK^X f{x),TxPKcYx) = df{pKS^))iTxPKcYx), 
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which shows that Xj{pK^{x)) = T^PKc^fix), because Uc is symplectic and T^prJ^x 
represents an arbitrary element of Tp^ (x)^- 

(b) To find / we choose a vector field X that is pi^^-related to Xj according to 
lemma 3.8 and then use 3.10, (b). ■ 

Definition 3.12. 

= {/ E I /|Me = 0}, 

= {/ G i?^^) I {f,h} € &\ for all h G i?^'^)}. 
Elements of I^^^ are called first class constraints. 

Lemma 3.13. (a) is a Poisson ideal of R^'^\ that is, it is an ideal of the ring 
such that if f e then {/, h} € /(^\ for all h G R^^h 

(h) Zj(c)R^^^ is a Poisson suhalgehra of R^^\ 

Proof, (a) Let f,g£ I^^^ and h S R^^\ Then it is immediate that f + g and kg 
belong to For any h G R^''\ we have {/, h]\Mc = Xh{f)\Mc = 0. 

(b) Follows from (a), using basic Poisson algebra arguments. ■ 

Lemma 3.14. The following conditions are equivalent for a function f G R^^\ 

(i) f e Zj,.,R(^\ 

(ii) f\Mc is locally constant. 

(Hi) Xf{x) G kera;c(x) for x G Mc- 

Proof. Assume (i). Then {/, h}\Mc = for ah h G -R^^), that is, df{x)Xh{x)\Mc = 
0. By lemma 3.10, (b), we know that Xh{x) represents any vector in T^Mc- We can 
conclude that f\Mc is locally constant, so (ii) holds. Now we will prove that (ii) 
implies (iii). Let f\Mc be locally constant. Then for all g G R^^'^ and all x G Mc, 

= Xg{f){x) = n{x){Xf,Xg){x) = ujc{x){Xf{x),Xg{x)). 

Since, again by lemma 3.10, Xg{x) represents any element of TxMc, we can conclude 
that Xf{x) G kerc<;c(a:), so (iii) holds true. Now we will prove that (iii) implies (i). 
Assume that Xf{x) G keTUJc{x), x G Mc- Then for all g G i?^*^^ and all x G Mc, 

{g, f}{x) = n{x){Xg,Xf){x) = u;c{x){Xg,Xf){x) = 0, 

that is, {g,f} G 1^"^^ Using this and the definitions, we see that / G Zj(c)R^''\ ■ 

Lemma 3.15. The map (pKc)* ■ Fi^^^ ~^ J~{Mc) defined in lemma 3.11 is a surjec- 
tive Poisson map and its kernel is therefore there is a natural isomorphism of 
Poisson algebras {pkJ^Kc) : R^^^^'' J='{Mc)- 

Proof. Surjectivity of (pKc)* ^tnd the fact that its kernel is I^^'^ follows immediately 
from lemma 3.11 and the definitions. This implies that {pKc)*ii':) is an algebra 
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isomorphism. Also, using the definitions, for f,g& R^"^^ and any x G Mc we can 
prove easily that 

{f,g}{x) = n{x){Xf{x),Xg{x)) = u:,{x){Xf{x),Xg{x)) 

= ^c{pkAx)) {Xf{pK,ix)),Xg {pkAx))) = {f,g}{pKAx)), 

where / = {pKc)*f^ 9 = {PKc)*9- Denote by -Kj^c) : R^'^^ R^'^'^/I^^^ the natural 
homomorphism of Poisson algebras. Then from the previous equalities we obtain 
{pkJ^i(c) {Trj{c){f),TTj(c){g)} = {f,g}, which shows that {pkJ^Hc) is a Poisson iso- 
morphism. In other words, we have the commutative diagram 

{PKc>* Nv I 

All the arrows are defined in a natural way and they are surjective Poisson algebra 
homomorphisms. ■ 

Equations of Motion and Physical Variables. It is immediate to see from the 
definitions that 

ker a;(x) n T^Mc ^ kercodx), (3-5) 

for all x € Mc- 

From now on we will assume the following. 

Assumption K2. (a) keru}{x) is a regular distribution, that is, it determines a 
regular foliation K and the natural projection px ■ M — )■ M, where M = AI/K, is 
a submersion. 

(b) The distribution kera;(3;) n TxMc is a distribution of constant rank. 

Theorem 3.16. The distribution kerct;(x) PI T^Mc is regular and has rank d^^\x). 
Its integral manifolds are SOMc, where S is an integral manifold o/kerw. Moreover, 
these integral manifolds give a foliation Kc of Mc which is regular, that is, the natural 
map pf.^ : Mc — > Mc, where Mc = Mc/Kc is a submersion. Besides, each leaf of the 
foliation Kc is foliated by leaves of Kc, which gives a naturally defined submersion 
Pj^ ^ : Mc Mc- In other words, we obtain the commutative diagram 

Mc-^Mc 




where each arrow is a naturally defined submersion. 
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Proof. The first assertion, about the rank of the distribution ker u}{x)r]TxMc, is easy 
to prove. Let xq € Mc- Then there exists a uniquely determined integral manifold 
S of the distribution kerw such that xq € S. Using that, by assumption, kera;(a;) fl 
TxMc is a distribution of constant dimension and that dim (kera;(x) n T^Mc) = 
dim {TxS n TxMc) we can conclude that the intersection S H Mc coincides with the 
integral leaf of the integrable distribution of 'ker ujr\TMc containing xq. So we obtain 
the foliation of Mc- Using (3.5) we can deduce that each leaf of the foliation K,. 
is foliated by leaves of Kc- The rest of the proof follows by standard arguments. ■ 

Lemma 3.17. (a) The following diagram is commutative 



M ^M- 



■TM 



Mr 



PKc 



fc 



■Me 



FcOtc 



— TM\Mr 




■ PcTM 




PKcKc 



where the arrows are defined as follows. The maps px, Vk^, Pj^^ cLnd Pj^^^c 
defined in Assumption Ki, Assumption K2 and theorem 3.16. By definition, the 
map fc is the inclusion. The map fc is an embedding defined by fc{S D Mc) = S, 
where S is a leaf of the foliation K. We will think of fc as being an inclusion. The 
vector bundle TM\Mc is the tangent bundle TM restricted to Mc. In other words, 
since fc is an inclusion, TM\Mc is identified via some isomorphism, called €c, with 
the pullback ofTM by fc. We call Fc the natural map associated to the pullback. 

(h) The presymplectic form lo on M passes to the quotient via px giving a 
uniquely defined symplectic form Co on M , satisfying pj^oo = oj. The presymplectic 
form ujc, which, by definition is fc^J, defines uniquely a presymplectic form Coc on 
Mc via pf^^ satisfying p*-. Cjc = ^c, = ft^- The energy function £ on M satis- 
fies d£{x)\kei[ uj{x) = 0, for all x € Mc, therefore it defines uniquely a 1-form on 
TM\Mc, called {Fc ° e)*d£ € F ( {TM\Mc)* I • Since £ is constant on each leaf of 



Kc, it also defines a function £c on Mc. Since TMc ^ TMc via the inclusion Tfc 
we have {Fc o e)*d£\TxMc = d£c{x), for all x € Mc. 

(c) Equation of motion (3.4) on Mc passes to the quotient Mc as 

cd{x){X{x), ) = {Fcoeyd£{x), (3.6) 

where X{x) € TxMc. This means that if X{x) G TxMc is a solution of (3.4) then 
X{x) := TxPj^ X{x), where x = Pj^ {x), is a solution of (3.6). Therefore, a solution 
curve x{t) of (3.4) projects to a solution curve x{t) = Pj^^{x{t)) of (3.6) on Mc. 
Equation (3.6) has unique solution X{x) for each x € Mc. This solution also satisfies 
the equation 

Cjc{x){X{x), ) = d£{x). (3.7) 

However solutions to equation (3.7) are not necessarily unique, since kerLDc(i) is 
not necessarily 0. 
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(d) The restriction of the energy function £\Mc satisfies 

d{£\Mc){x)\keru}c{x) = 0, 

for all X € Mc, therefore there is a uniquely defined function £c on Mc such that 
Pkc^c = £\Mc. The equation 

u>,ix)iX{x), ) = d£,{x) (3.8) 

has unique solution X{x) for x € Mc- If X{x) is a solution of (3.6) then X[x) = 
TxPj^^f^^X{x) is a solution of (3.8). Therefore, a solution curve x{t) of (3.6) projects 
to a solution curve x{t) = pj^ ^ (^(^)) of (3.8) on Mc. 

Proof, (a) The equality px^ = Pj^ ^ ° Pk proven in theorem 3.16. The 

equahty Tj^joF^oe^ = fc°Tff results immediately from the definitions. The equality 
Pk ° fc = fc ° Px results by applying the definitions and showing that, for given 
X G Mc, Pk ° fc{x) = Sx = fc °Pk (^)) where Sx is the only leaf of K containing x. 

(b) Existence and uniqueness of a) and is a direct consequence of the definitions 
and standard arguments on passing to quotients. For any x € Mc we know that 
there exists a solution X of equation (3.4), from which it follows immediately that 
d£{x) I ker oj[x) = 0. The rest of the proof of this item consists of standard arguments 
on passing to quotients. 

(c) We shall omit the proof of this item which is a direct consequence of the 
definitions and standard arguments on passing to quotients. 

(d) If X{x) is a solution of uj{x){X{x), ) = d£{x) then it is clear that it also 
satisfies uJc{x){X{x), ) = d{£\Mc){x). It follows that 

d{£\Mc){x)\keTU}c{x) = 0, 

for all X G Mc- The rest of the proof is a consequence of standard arguments on 
passing to quotients. ■ 

Remark. Recall that the locally constant function d^'^\x) on Mc is the dimen- 
sion of the distribution kera;(3;) n T^Mc on Mc and also the dimension of the fiber 
of the bundle S'^''\x). If d^"'^ (x) is nonzero then there is no uniqueness of solution to 
equation of motion (3.4), since solution curves to that equation satisfy, by definition, 

u}{x){±, ) = d£{x)\TxM, 

where (x, x) G TxMc. Passing to the quotient manifold Mc eliminates this inde- 
terminacy and uniqueness of solution is recovered. This is related to the notion of 
physical variables mentioned by Dirac. 

3.3 First Class and Second Class Constraints and Constraint Sub- 
manifolds. The Tangent Bundle to a Second Class Subman- 
ifold. 

As we have said before, an important topic in the theory of constraints as devel- 
oped by Dirac is the distinction between first class and second class constraints. 
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His treatment is intended to solve systems with constraints coming from degenera- 
cies in the Lagrangian from both the classical and the quantum mechanics point of 
view. The Poisson algebra structure of functions on a symplectic manifold is the 
context in which this theory is developed and it is not very geometric and almost 
no attention is paid to the constraint submanifolds defined by the several equations 
involved. Among several interesting references we cite Sniatycki [1974] which has 
several points of contact with our work. 

In this paragraph we shall give a geometric framework and describe its close 
relationship to the Poisson algebra point of view to deal with the notions of first 
class and second class constraints and functions and also first class and second 
class submanifolds. These notions only depend on the final constraint submanifold 
Mc and the ambient symplectic manifold and do not depend on the primary 
constraint Mq = M or the Hamiltonian H : N ^ M. Accordingly, in this paragraph 
we will adopt an abstract setting, where only an ambient symplectic manifold and a 
submanifold are given. This kind of abstract setting was studied in Sniatycki [1974], 
in particular the notion of second class constraint submanifold and its connection 
with the Dirac bracket. 

Then we will go back to equations of motion in the next subsection, where the 
role of both the final and the primary constraint is essential. 

The definitions given at the beginning of this section inspire the following one. 

Definition 3.18. Let {P,fl) be a symplectic manifold and S Q P a given submani- 
fold. Then, by definition, 

RiS,P) .- {f e j-(p) I Xf{x) G T^S, for all x G 5} 
/(^'^) := {/Gi?(^'^) 1/15 = 0} 

Elements of R^^'^^ are called first class functions. Elements of I^^'^^ are called 
first class constraints. The submanifold S is called a first class constraint 

submanifold if for all f G J^{P) the condition /|5 = implies f G that is, 

I(S,P) — where I[s,P) 'is the ideal of the ring J-{P) of all functions vanishing 

on S. 

Obviously, using the notation introduced before in the present section, R^^'^'^^ = 
and iC^^-^^) = li"). All the properties proven for R^'^^ and I^'^^ hold in general 
for ii(^'-f') and For instance, T^S is the set of all Xf{x), f G R'^^'^\ Ev- 

ery function / G iJ^"^'^) satisfies df{x){Xr,) = 0, for all G ker{il.{x)\T^S) and 
keT{n{x)\T^S) is the set of all Xf{x), f G 

The following lemma 3.19 is one of our main results. It studies the vector 
subbundles V C TP\S which classify all second class submanifolds containing 
S" at a linear level, that is, is tangent to the second class submanifold. For such 
a second class submanifold, which is a symplectic submanifold, the Dirac bracket 
of two functions F and G at points x G 5 can be calculated, by definition, as 
the canonical bracket of the restrictions of F and G. This has a global character. 



3.3 First Class and Second Class Constraints and Constraint Submanifolds... 23 

A careful study of the global existence of a bracket defined on sufficiently small 
open sets U ^ P containing S which coincides with the previous one on the second 
class submanifold will not be considered in this paper. However, to write global 
equations of motion on the final constraint submanifold one only needs to know the 
vector bundle V^, which carries a natural fiberwise symplectic form. 

Local descriptions of second class constraint submanifolds on sufficiently small 
open sets give the usual expression of the Dirac bracket. 

All these are fundamental properties of second class constraints and constraint 
submanifolds, and theorem 3.20 collects the essential part of them. 

Lemma 3.19. Let (P, he a symplectic manifold of dimension 2n and S Q P a 
given submanifold of codimension r. Letuj he the pullhack ofQ to S and assume that 
keruj{x) has constant dimension. Assume that S is defined regularly hy equations 
<j)i = 0, . . . ,(j)r = on a neighhorhood U ^ S and assume that we can choose a 
subset {xi, ■■■ ,X2s} ^ {0i5---50r} such that det{{xi,Xj}{x)) "j^ for all x G S, 
where we assume that 2s = rank({0j, for all x £ S. We shall often denote 

'^^ji^) — {XiiXj}{x) iiT-d Cx{x) the inverse of cfj{x). Moreover, we will assume that 
the following stronger condition holds, for simplicity. Equations <j)i = Bi, . . . ,(f)r = 
Br and xi = Ci, . . . , X2s = define submanifolds of U regularly, for small enough 
Bi, . . . , Br and Ci , . . . , C2s ■ 
Then 

(a) 2s = r — dimkerw = 2n — dimS — dimkertj. There are V'fc G k = 
1, . . . ,r — 2s, which in particular implies {tpi^, ipi}{x) = 0, {ipk, Xi}{x) = 0, for k, I = 
1, . . . ,r — 2s, i = 1, ... ,2s, such that d'ipi{x), . . . , dtpr-2six) , dxi{x) , . . . , dx2s{x), are 
linearly independent for all x £ S. Moreover, X^j(x), . . . ,X^^_2^{x) from a basis of 
kevuj{x), for all x £ S and dipi{x), . . . ,dipr-2s{x),dxi{x), . . . ,dx2s{x) form a basis 

ofiT^sy. 

(b) The vector subbundle C TP\S with base S and fiber 

= span(X^,(x), . . . C T.P, 

satisfies 

V^r^T,S = {Q} 
V^®kfixuj{x) = (T^S) 
{V,^fn{keru;{x)f=T,S, 

xeS. 

(c) There is a neighhorhood U of S such that the equations xi = 0, . . . , X2s = 0, 
on U define a symplectic submanifold such that S Q S-^ and 

T,S^ = {V^f 
T,s^ e = T^P, 

for X £ S^, where we have extended the definition of Vx for x £ using the 
expression 

= span {X^^ (x), . . . , X^^^ {x)) C T^P, 



(3.9) 
(3.10) 
(3.11) 
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forx £ S^. The submanifold has the property I(^s,sx) ^ that is, S is a first 

class constraint submanifold of S^, defined regularly by ijji\S^ = 0, i = 1, . . . , r — 2s, 
and ipi\S^ G jis,sx) ^ ^ _ , , , — 2s. Moreover, it has the only possible dimension, 
which is dim 5-^ = dimS' + dimkerw = 2n — 2s, for symplectic submanifolds having 
that property. It is also a minimal object in the set of all symplectic submanifolds 
Pi ^ P, ordered by inclusion, satisfying S C Pi. 

(d) Let V be any vector subbundle ofTP\S such that 

y ekertj = (r5)^, (3.12) 

or equivalently, 

@{keTUjf = {TSy. (3.13) 

Then dim 14 = 2s, for x G S. Let be a submanifold of P such that T^S^ = V^, 
for each x € 5. Then S is a submanifold of . Such a submanifold always 
exists. Moreover, for such a submanifold there is an open set U ^ P containing S 
such that nU is a symplectic submanifold of P. 

Let X £ S and let x'l = 0, . . . , X2s = 6e equations defining PI U' for some 
open neighborhood U' P and satisfying that dxi{x),... ,dx2six) are linearly in- 
dependent for X e n [/'. Then, dxi{x), . . . , dx2six), dipi(x), . . . , d'il)r-2s{x) are 
linearly independent and de[,{{x'i,x'j}{x)) ^ 0, for x e fl U' . All the properties 
established in (a), (b), (c) for xi, - ■ ■ ,X2s on S hold in an entirely similar way for 
x'l, ... , X2s' on S nU' . In particular, nU' = . 

(e) Let bS^ be the pullback of Q to and {,}^ the corresponding bracket. For 
given F,G £ ^{P) define F^ := F — XiCxiXj^P} o,nd also 

{F, := {F, G} - {F, x.}c^^{x„ G}, 

which is the famous bracket introduced by Dirac, called Dirac bracket, and it is 
defined for x in the neighborhood U where cfj{x) has an inverse c^{x). Then, for 
any x G S^, 

{F^,Xi}{x)=0 

for i = 1, . . . , 2s, and also 

{F^,GJ(x) = {F,G}(^)(x) =o.^(x) {Xf\s4x),Xg\sAx)) = {F\S\G\S^Ux). 

If we denote l^he Hamiltonian vector field associated to the function F £ 

J-{P), with respect to the Dirac bracket {, then the previous equalities are equiv- 
alent to 

^[x)A^) = ^F^ix) = Xpisx{x). 
The Jacobi identity is satisfied for the Dirac bracket {F, on S^, that is, 

{{F, , (x) + {{H, F}(^) , (x) + {{G, , F}(^) (x) = 0, 

for X £ S^. 

(f) Let U be an open neighborhood of S such that c^ (x) is invertible for x £U. 
For each C = (Ci, . . . , C2s) G K^** let x't = Xi~ Ci and define S^^ by the equations 
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Xi{x) = 0, i = l,...,2s, X £ U. For any C in a sufficiently small neighborhood 

c c 
of 0, is a nonempty symplectic submanifold of P. Define the matrix cfj (x) = 

{xf\xf}{x)! (^^d also c^c{x) as being its inverse, x € f7. Then, the equalities 

cfjix) = {x^,XJ}{x) = {X?,xf}{x) = cf/(x), (3.14) 

and also, 

{F,G}(^)(x) = {F,G}(^c)(x) (3.15) 

are satisfied for all x (z U . All the definitions and properties proved in (e) for the 
case C = hold in general for any C in a neighborhood of small enough to ensure 

(J 

that is nonempty. In particular, the equalities 

{F^c,G.^c}{x) = {F,G}(^^c){x) = w^''(x) {Xp^g^c{x),X^^g^c{x)^ 

= {F\S^'' ,G\S^''}^c{x) (3.16) 

and 

hold for X G , and any C in such a neighborhood. The Dirac bracket {F, 
satisfies the Jacobi identity for F,Gg J~{U) and the symplectic submanifolds 
are the symplectic leaves of the Poisson manifold {J-{U),{, }{x))- By shrinking, if 
necessary, the open set U and for C in a sufficiently small neighborhood ofOG M^*, 
the equations ipklS'^ =0, k = I, . . . ,r — 2s, define regularly a first class constraint 

submanifold S'-" C S-^ '^^ U, and the functions V'fcl'S'^ € R^^ '^^ C F[S^ ) are 
first class constraints, that is, ipk\S-^'^ G jis'-'\sx ) = . . . ^7. _ 2s. We have that 
dimS'^ = dimS''^ + dimkero;'^, where uj'^ is the pullback of to S'^' . One has 
dimS''-^ = dimS and dimS^ = dimS"^, therefore dimkerc<; = dim ker c<;'-" . 

Proof, (a) Let x S. We are going to use lemmas and corollaries 3.1-3.7 with 
E := TxP; V := T^S] 7^ := d(j)i{x), i = 1, . . . ,r] Q := Q{x); u := oj{x); f3 = 0. 

Elements i/j = ^'(pi, A* € J'(P) such that ■0 S which implies X^{x) £ T^S 

for X £ S, must satisfy {ip,(j)j}{x) = 0, or, equivalently, X^{x)dcj}i{x) (^X^_^{x)) = 
0, for j = l,...,r, X S 5. Using lemma 3.6 we see that X^{x) E kera;(x). 
Since one can choose r — 2s linearly independent solutions, say Aj = {Xj, . . . , A[), 
i = 1, ... ,r — 2s, we obtain elements V'j S namely, 0j = Xl(j)j, such that 

((iV'i(x), . . . , (iVr-2s(a;)) are linearly independent, or, equivalently, taking into ac- 
count lemma 3.7, that (^^^(x), . . . , X^^_^^{x)) is a basis of kera;(x) for x £ S. If 
dipi{x), d'ipr-2six), dxi{x), dx2six), were not linearly independent, then 
there would be a linear combination, say x = with at least one nonzero coeffi- 
cient, and some x £ S, such that ^^(2;) = fi^dipk{x) for some k = 1, . . . ,r — 2s. 
But then, for any j = l,...,2s, {x,Xj}{x) = dx{x)Xx.{x) = fi''dipk{x)Xx^{x) = 
—^^dxj{x)X^^{x) = 0, which contradicts the fact that det {{Xi^Xj}{x)) 0. Using 
this and the fact that V'i = 0, = 0, i = r — 2s, j = 1, . . . , 2s define S regularly, we 
can conclude that dipi{x), . . . , d'4jr-2s{x), dxi{x), . . . , dx2s{x) is a basis of (TxS)° . 
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(b) If X'X^^{x) e T^S then dxj{x)\' X.^^{x) = 0, j = 1, . . . , 2s, which imphes 
X^{Xj-:Xi] = 0, j = 1, . . . , 2s, then A* = 0, i = 1, . . . , 2s, which proves (3.9). To 
prove (3.10) we apply the operator ^ to both sides and obtain the equivalent equality 
span((ixi(a^), • • • , dx2s(2:)) © span((i'0i(a;), . . . , (iV'r-2s(a^)) = {TxS)°, which we know 
is true, as proven in (a). To prove (3.11) we apply the orthogonal operator ^ to 
both sides of (3.10). 

(c) Since dxi{x), . . . ,dx2s{x) are linearly independent for x G 5 they are also 
linearly independent for x in a certain neighborhood U of S. Then Xi(a;) = 
0, . . . ,X2s{x) = define a submanifold of U containing S. To see that it is 
a symplectic submanifold choose x (z and apply corollary 3.7 with E := T^P; 
V := TxS^; j3 := 0; 7^ := dxi, i = l,...,2s; oj := Q,{x)\TxS^. We can con- 
clude that dim(ker r2(x)|Ta;S'^) = 0. Now let us prove that TxS^ = {V^y\ namely, 
TxS>^ = span(dxi(x), . . .,dx2s{x))° = {{V^ff = {V^f. FVom this, using that 

is symplectic one obtains T^S^ © Vx = TxP. To prove that 5" C 5-^ is a first class 
constraint submanifold defined by first class constraints ipi\S-^, i = 1, . . . , r — 2s, on 
S^, we observe first that it is immediate that ipi\S^ = 0, i = 1, . . . , r — 2s, define S 
regularly. It remains to show that X^.\s^{x) G TxS, i = 1, . . . ,r — 2s, x G S , where 
X^.^gx is the Hamiltonian vector field associated to the function tpilS^ with respect 
to the symplectic form u}^. This is equivalent to showing that 

X^^\sx{x) {ijj\S^) = 0, 

for X G 5 or, equivalently, 

LV^ix) (^X^,^|5x(x),X^.|5x(x)) = 0, 

for X G 5. We know that is the pullback of to 5^ and ipilS^ is the pullback of 
ipi, via the inclusion 3^(^11, then we have 

uj^ix) (x^^|5x(x),X^^.|5x(x)) = n{x) (X^,(x),X^^,(x)) = 0, 

for X G 5", since ipi are first class constraints, i = l,...,r — 2s. Finally, using 
the definitions we can easily see that dim = 2n — 2s and that dim S = 2n — r 
and from (a) we know that dimkerw = r — 2s. We can conclude that 2n — 2s = 
dim 5 + dim ker u . 

(d) We know that, for x G S", dimT^-S = 2n — r, and dimkera;(x) = r — 2s; 
then using (3.12) we obtain dim 14 = 2s. Also from (3.12) we immediately deduce 
applying ^ to both sides, 

n {kevujf = TS, 

in particular TS C V^. Let g he a given Riemannian metric on P and let Wx be 
the (7-orthogonal complement of TxS in V^, in particular, Wx ®TxS = V^, for each 
X G S. Define 

= {exp{twx) I Wx G Wx,gix){wx,Wx) = 1, |t| < r(x),x G S}. 

By choosing t(x) appropriately one can ensure that is a submanifold and, more- 
over, it is easy to see from the definition of 5^ that TxS^ = Wx (BTxS = V^, for 
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each X £ S. We leave for later the proof that PI C/ is a symplectic submanifold 
of P, for an appropriate choice of the open set U, which amounts to choosing t{x) 
appropriately. 

Assume that dx'i{x) , . . . , dx'2s{x) are linearly independent for x G n U' . 
Since {dXi{x),V^) = {d'x'^{x),TxS^) = for x S 5 and i = l,...,2s, we can 
deduce that d'x'^{x) G {Vx) i that is, dx'i{x) € V^. Then using (3.13), we see 
that dx'i{x), . . . , dx'2s{^)-> d'4'i{x)-, • • • , dipr-2six) are linearly independent and span 
V^Qikevuf = {TSy. If det({x-, XjK^;)) = for some x € 5 then A^x-, XjK^;) = 0, 
where at least some A* / 0, i = l,...,2s. Let X^x'i = "^i then {99, = 0, 

J = l,...,2s. On the other hand, since ip\S = 0, then {Lp,il;j}{x) = 0, j = 
l,...,r — 2s. We can conclude that Lp € /(■^''^ ) and then Xip{x) € kera;(x), in 
particular, = n^X^.{x), which implies X^dXiix) = dip{x) = fi^diljj{x), contra- 

dicting the linear independence of dxi{x), . . . , dx'2s{x), dipi{x), . . . , dipr-2s{x). 

It follows from which precedes that by replacing Xj by x^, i = 1, . . . , 2s and S 
hy S CiU' all the properties stated in (a), (b) and (c) are satisfied. In particular, 
S nU' = and S OU' is symplectic. It is now clear that by covering S with open 
subsets like the U' we can define U as being the union of all such open subsets and 
one obtains that S DU is a symplectic submanifold. 

(e) Let X G S^. Then, since F-^ = F — Xi^x {Xj ■• ■• obtain 

{Fx,Xk]{x) = {F,Xk]{x)-{x^,Xk}{x)4{x){xj,F}{x) = {F,Xk}{x)+{Xk,F}{x) = Q. 
Using this we obtain 

{F^,GJ(x) = {F^,G-Xfc4'{x/,G'}}(x) = {F^,G}(x) 

= {F,G}(x) - {x.,GKi{x„i^} = {F,G\x){x). 

For any F e ^{P), x G and A; = 1, . . . , 2s, we have Xp^ ix)xk = {Xfc; -^x}(^) ~ 
0, so Xp^{x) G TxS^. Therefore, for any G T^S^, 

uj^ix) {Xf^{x),Yx) = n{x) {Xf^{x),Yx) = dF^ix)Yx = d{F^\S'>') (x)y, 

= d{F\S^) {x)Yx = io^x) {Xp^sx{x),Yx) , 

which shows that Xf^{x) = Xp^^xix), where both Hamiltonian vector fields are 
calculated with the symplectic form u^. Using this, for any G € J^{P) and any 
X G S^, one obtains 

{G^,Fj(x) = Xf^{x)G^ = Xp^{x)G = Xf\sx{x)G\S^ = {G\S^ ,F\S^}^{x). 

The equality X(^^'^^p{x) = Xp^{x) = Xp\sx{x) is an immediate consequence of the 
previous ones. The Jacobi identity for the bracket { , }(^) follows using the previous 
formulas, namely, for x S S"-^, one obtains 

{{F,G}(^),F}(^)(.x) = {{F,G}^^)\S^,H\S^Ux) = {{F\S\G\S^}^, H\S^Ux), 

where the bracket in the last term is the canonical bracket on the symplectic manifold 
S^, for which the Jacobi identity is well known to be satisfied. 
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(f) The equalities (3.14) and (3.15) are proven in a straightforward way. The 
equations (3.16) and (3.17) follow easily using a technique similar to the one used 
in (e). Using all this, the proof of the Jacobi identity for the bracket { , }(^) on U 

goes as follows. Let x & U and let C be such that x E 5-^ . For F,G,H G J^{U) 

c 

using (e) we know that the Jacobi identity holds for { , }(;^C) on S-^ . But then, 
according to (3.15) it also holds for { , }(^) for all x € 5^ . Now we will prove 
that S^^ are the symplectic leaves. Since they are defined by equations xf = Oi 
i = I, ... ,2s on U we need to prove that {F, xF}x(^) = 0, x G S'~^ , for all F G J^{U), 
i = 1, ... ,2s. Using (3.15) and (3.16) we see that {i^, X?}(x) (^) = {F,X?}ix'^){x) = 
{F\S^ ,Xi\S^ }xc(^) ~ 0. To finish the proof, observe first that, since Xi = 0) 
V'i = 0, i = 1, . . . , 2s, j = 1, . . . , r — 2s define regularly the submanifold 5 C [/, by 
shrinking U if necessary and for all C sufficiently small, we have that xf = 0, = 0, 
i = 1, . . . ,2s, j = 1, . . . , r — 2s define regularly a submanifold C U and therefore 
il^j\S'^ = 0, j = 1, . . . , r — 2s, define regularly S as a submanifold of . To prove 
that it is a first class constraint submanifold and that 'ilJj\S'^ , j = 1, . . . ,r — 2s, 

are first class constraints, that is tl^j\S^^ G \ j = 1, . . . ,r — 2s, we proceed 

in a similar fashion as we did in (c), replacing X by X • The fact that ipjl^"^ ; 
j = l,...,r — 2s, are first class constraints defining 5*^ implies that dimS^^ = 
dim 5*^ + dimkerw*-". From the definitions one can deduce that dimS"^ = dim 5 
and dim S*-^ = dim S"-^, therefore dim kertj = dim ker a; . ■ 

The following theorem summarizes the essential part of the previous lemma. 

Theorem 3.20. Let {P, Q) be a symplectic manifold, S C P and let u be the pullback 
of Vt to S. Assume that kercj has constant rank. Let V be a vector subbundle of 
TP\S such that F kercj = {TS)^. Then there is a symplectic submanifold 
containing S of dimension dimS + dimker a; and the condition T^S^ = V^, for all 
X G S" holds. The vector bundle is called the second class subbundle tangent 
to the second class submanifold . For given functions F, G on P, one defines the 
Poisson bracket {F,G}^(x) := {F\S^ ,G\S^}{x), x G 5. We call this the V-Dirac 
bracket on S. 

Remark (a) Since we are interested mainly in describing equations of motion we 
will not consider the definition of a global Poison bracket on a neighborhood of S 
such that one of its symplectic leaves coincides with . By choosing a local regular 
description Xfe = of one obtains the usual expression for the Dirac bracket, as 
it was shown in lemma 3.19. Under our strong regularity conditions the symplectic 
leaves of the Dirac bracket give a (local) regular foliation of a neighborhood of the 
final constraint submanifold S. This implies by the Weinstein splitting theorem 
(Weinstein [1983]) that there are local charts where the Dirac bracket is constant. 

(b) The tangent second class subbundle in a sense (modulo tangency) clas- 
sifies all the possible second class constraint submanifolds containing a given sub- 
manifold S C. P. It carries enough information to write the Dirac brackets along the 
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final constraint submanifold S and therefore also equations of motion, as we show 
in subsection 3.4. 

3.4 Equations of motion 

We are going to describe equations of motion in the abstract setting of the para- 
graph First Class and Second Class Constraints and Constraint submanifolds of 
subsection 3.2, that is, a symplectic manifold {P, and a submanifold /S C P, 
defined regularly by equations (pi = 0, i = 1, . . . ,a. We are going to assume all the 
results, notation and regularity conditions of that paragraph. We need to introduce 
in this abstract setting, by definition, the notions of primary constraints, primary 
constraint submanifold and the Hamiltonian. 

The primary constraint submanifold is a given submanifold S' Q P con- 
taining S, and in this context, S will be called the final constraint. We will assume 
without loss of generality that S' is defined regularly by the equations (pi = 0, 
i = 1, . . . , a', with a' < a, where each (pi, i = 1, . . . ,a' will be called a primary 
constraint while each (pi, i = a' + 1, . . . ,a will be called a secondary constraint, 
for obvious reasons. In this abstract setting the Hamiltonian is by definition a 
given function H G J~iP)- 

The equations of motion can be written in the Gotay-Nester form, 

Q{x){x,Sx) = dH{x){Sx), 
where {x,x) G T^S, for all 6x € T^S' . 

Now we will transform this equation into an equivalent Poisson equation using 
the Dirac bracket. 

The condition {H,^p}(x) = 0, for all x G 5 and all first class constraints ^p will 
appear later as a necessary condition for existence of solutions for any given initial 
condition in S, so we will assume it from now on. 

The total Hamiltonian is defined by Ht = H + X^(pi, i = 1, . . . ,a' where the 
functions A* G C°°{P), i = I,..., a must satisfy, by definition, {Ht, (pj}{x) = 0, 
_7 = 1, . . . , a, X G 5 or, equivalently, 

{H,(pj}{x) + X'{(pi,(pj}{x) = 0, x e S,i = l,...,a', j = l,...,a, 

sum over i = 1, . . . ,a' . We assume that the solutions (A"*^, . . . , A"') form a nonempty 
affine bundle A — > 5". 

For each section of A one has a Hamiltonian Ht{x) = H{x) + X'{x)(pi{x), x £ P, 
and an equation of motion 

= {dHrf. 

Equations of motion can be described nicely using the Dirac bracket as we will 
see in a moment. Choose first class and second class constraints 

(V'l,...,V'a-2s,Xl,---,X2s) 
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as in lemma 3.19, then since {-ffT)Xi}(^) = 0) ^or any function F € J~{P) and any 
X G S we obtain 

{Ht,F}^^){x) = {Ht,F}{x) - {HT,Xi}c^{x){Xj,F}{x) = {Ht,F}{x), x e S. 

Then we can describe equations of motion in terms of the Dirac bracket as 

XhA^) = X^^)^Hri^), xeS. (3.18) 

We want a more precise description of equations of motion. The total Hamiltonian 
has the equivalent expression 

HT = H + X"4 + fi'^x'j, (3.19) 

where ip'^ G I{S',P) ^ R^^'^\ i = 1, . . . ,a' — s' are such that {dip[{x), . . . , dip'^,_^,{x)) 
form a basis of 

{#(x) I V€/(s',P)ni?(^'^)}, 

for all x £ S, while x'i ^ ^{S',P)^ i = ^, ■ ■ ■ , s' , are such that 

(dil^'^ix), . . . ,d'il)'^,_^,{x),dxi{x), . . .,dx's'{x)) 

form a basis of 

{dil^ix) I ip € I{s',P)} 

for all X G S. Then {il>[, . . . , V'l'-s'i x'l) • • • i x's') can be chosen as a set of primary 
constraints which justifies the expression (3.19) for the total Hamiltonian. One can 
see that the rank of the matrix {xj ,x'i}i^)^ j = 1, . . . , 2s, i = 1, . . . , s', x € is 
s' . Now the conditions (3.18) are equivalent to {HT-,'4'j}[x) = 0, j = 1, . . . , a — 2s, 
which gives {H, ijjj}{x) = 0, j = 1, . . . , a — 2s and {Ht, Xi}i^) = 0, i = 1, . . . , 2s, 
for all x G S* which gives 

{H,Xi}{x) + fi'^{x){x',,Xi}{x) = 0, 



i = 1, . . . , 2s, X E S, from which we obtain ^'^ , j = 1, . . . , s', as a well-defined 
function (on a neighborhood of 5, then we can extend it arbitrarily to P). One can 
write Xj = V'fc + b^jXi with a^, h^- G F{P) uniquely defined on 5, for j = 1, . . . , s', 
A: = 1, . . . , a — 2s, / = 1, . . . , 2s. Then the total Hamiltonian can be written 



Ht = H + A'Vi + li''a]^l^k + f^''b'jxi 
= H + iP^s,S',H) + A'V^ + ^^'^b^jXl, 

where ip(s,S',H) = fJ-''' CLji^k is a first class constraint, '4'{S,S',H) £ I^^''^^- 

We can conclude, using the fact that {xfc, -^}(x) (^) = 0> for all x G 5 and 
A: = 1, ... ,2s and (3.18), that 

{Ht,F}{x) = {Ht,F}^^^{x) 

= {i7 + V'(s,5',i/) + A'V:,F}(^)(x) 
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for any function F € J~'{P) and any rr € 5*. 
We shall cah 

H{x,S,S') ■■= H + tp(s,S',H) + >^"i^[s,S',i) 
the Hamiltonian of the system with respect to the Dirac bracket { , }(^x)^ where we 
have denoted ^p'^ by tp'f^g gi j) to emphasize that these functions depend on S and S'. 
We have proven the following theorem. 

Theorem 3.21. Let {P,0,) be a symplectic manifold and S Q S' Q P given regu- 
larly defined submanifolds. In the situation described above, equations of motion on 
S can be written in the following equivalent ways: 

(a) 

Xhj, (x) = X(^^)^H^^^s,s') = ^ix),H+ij(s,.s',H)+y'^- (3-20) 
for all X G S. 

Here H € J~{P) satisfies, by definition, {H,'ip}{x) = 0, for all x G S and all 
primary first class constraints ip; ip(s,S',H) a first class constraint depending on 
S, S' and H; {V'l, • • • , "^a'-s'l ^ hs',P) ^ R^^"^^ is a maximal independent (that is 
{dip'i{x), . . . ,dip'^,_^,{x)} is linearly independent for each x G S) set of first class 
primary constraints and A'*, i = 1, . . . ,a' — s' are arbitrary parameters. There is 
uniqueness of solution if and only if there are no first class primary constraints, that 
is, a' — s' = 0. 

(^) 

F = {F,H(^^^s,s')}{x)- 



Q{x){x,6x) = dH{x){5x), 
where {x,x) € T^S, for all 5x G T^S' . 

Remark, (a) The family of Hamiltonian vector fields (3.20), each one of them 
tangent to S, depends on a finite number of arbitrary parameters A'* G M, z = 
1, . . . , a' — s', and it generates an affine distribution on S. This should be compared 
with the affine bundle in equation (2.7). Any vector field, even time-dependent, Xt 
on P whose restriction to /S is a section of that distribution gives equations of motion. 
Note that, since il^i\S = 0, then for any choice of functions, even time-dependent, 
\'l e J^{P), i = 1, . . . ,a' — s' we have 

K'^{x),^'S^) = ^(x),xf^'S^)^ 

for all x G S". 

(b) The equations of motion can be globalized, using the bracket { , }^, as far 
as one can find a global '>p(s,S',H)- 
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4 Future work 

The present work should be followed immediately by an extension of the Dirac 
theory, and also the dual Gotay-Nester theory for the case of a Dirac dynamical 
system (x, x) © d£{x) G D^. This will expand the field of applications, for instance, 
one will have a theory of constraints for LC circuits and holonomic systems, if the 
Dirac structure D is integrable. 

The reduction theory for the constraint algorithm will be also the purpose of 
future work. This can be approached using the category of Dirac anchored vector 
bundles. This will extend part of the results in Cendra, Marsden, and Ratiu [2001]. 

Singular cases, where the strong regularity assumptions made in the present 
paper are weakened in several ways are also important and will be the purpose of 
future work. 



A Appendix 

Lemma A.l. There is a canonical inclusion (p: TQ © T*Q — )■ T*TQ. In addition, 
consider the canonical two-forms ujt*q and ujt*tq on T*Q and T*TQ respectively, 
the canonical projection P'Ct*q '■ TQ@T*Q T*Q, and define the presymplectic two- 
form uj = pr^,Q ojt*q on TQ(BT*Q. Then the inclusion preserves the corresponding 
two-forms, that is, u = ip*ujT*TQ- 



Proof. Utq:TQ^Q and ttq : TTQ 
consider the dual tangent rhombic 



TQ are the tangent projections, we can 



TTQ 




Define (p:TQ®T*Q 



T*TQ by ifivg © ag) e T*TQ, 
Lp{vq © aq) ■ = ag ■ Ttq{w^^), 



for Wi,^ G Ty^TQ. Here Vq © aq denotes an element in the Pontryagin bundle over 
the point q (z Q. Note that the following diagram commutes. 



TQ®T*Q- 

P^TQ 

TQ 



Id 



TQ 



T*TQ 

TTTQ 

^TQ 



Let us see that is an injective vector bundle map from the bundle pr^^g : TQ © 
T*Q — )■ TQ to the cotangent bundle tttq- T*TQ — )■ TQ, over the identity of TQ. 
The last part of this assertion follows from the commutative diagram above. 
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First, if ^p{vq © a^) = ^p{v'^, © a'^,) then both sides are in the same fiber T*^TQ = 
T*, TQ, so Vq = v' ,. Also, for all w^^ G T* TQ we have 

(p{Vq © aq) ■ Wy^ = ip{Vq © Q^) • Wy^ 

or 

aq ■ Ttq{w^^) = a'g ■ Ttq{w^^). 

Since tq : TQ — > Q is a submersion, we have aq = a'q and (p is injective. 
Second, is linear on each fiber, since 

^{vq®{aq + \a'q))-w^^ = {aq + Xa'^) ■Ttq{w^^) = ip{Vq®aq)-'Wi,^ + Xip{Vq®a'q)-W^^ 

For the second part of the lemma, let us recall the definition of the canonical 
one-form on 9t*q € Q^{T*Q). For aq G T*Q, 6T*Q{aq) is an element of T*^T*Q 
such that for £ TagT*Q, 

Or-Qiaq) ■ Wag = aq{TTTQ{Wag)), 

where ttq : T*Q — )■ Q is the cotangent bundle projection. 

With a similar notation, the canonical one- form 9t*tq G 0,^ {T*TQ) is given by 

Or-TQioivg) ■ Wa,^ = a^g{TTiTQ{Wa^^)). 

Pulling back these forms to the Pontryagin bundle by ip and the projection pr^^.g : TQ® 
T*Q T*Q, we obtain the same one-form. Indeed, for w G T^^^a^iTQ © T*Q), we 
get on one hand 

{pr^*QOT*Q){Vq eaq)-W = eT*Q{aq) ■ TpTt,q{w) = aq ■ T{tTq O ^Ij„q){w), 

and on the other hand 

{ip*0T*TQ){Vq ®aq)-W = OT*TQ{^{Vq © aq)) ■ Tip{w) = 
ip{Vq © aq) ■ T-KTQ(Tip{w)) = aq ■ T{tq O TTtQ O 

However, the following diagram commutes 

TQ © T*Q ^ T*TQ 




so TTQ o picj^tQ = Tq o tttq o ip and therefore pr^^tg Ot*q = (p*9t*tq- Taking (minus) 
the differential of this identity, we obtain lo = ip*ujT*TQ- ^ 

In local coordinates, if we denote the elements oiTQ®T*Q and T*TQ by (g, v,p) 
and {q,v,p,v) respectively, then ip{q,v,p) = {q,v,p,0). 
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